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We present a theoretical framework for evaluating effective interactions between localized spins 
mediated by itinerant electrons in double-exchange models. Performing the expansion with respect 
to the spin-dependent part of the electron hopping terms, we show a systematic way of constructing 
the effective spin model in the large Hund’s coupling limit. As a benchmark, we examine the 
accuracy of this method by comparing the results with the numerical solutions for the spin-ice type 
model on a pyrochlore lattice. We also discuss an extension of the method to the double-exchange 
models with Heisenberg and XY localized spins. 
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I. INTRODUCTION 


Spin-charge coupled systems, which consist of itiner¬ 
ant electrons interacting with localized moments, are of 
special interest in the study of itinerant magnets for 
their rich physics from the interplay between magnetic 
moments and itinerant electrons. A key aspect that 
gives rise to the rich physics is the effective interac¬ 
tions between localized moments mediated by the itin¬ 
erant electrons. In the weak coupling limit, the itinerant 
electrons induce long-ranged effective exchange interac¬ 
tions with oscillating signsA— The weak coupling theory 
of effective spin interactions, so-called Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interactions, has achieved a suc¬ 
cess in the study of metallic magnets, such as transition- 
metal^ and rare-earth^ compounds and spin-glass behav¬ 
ior in alloys lightly doped with magnetic ionsi^ 

On the other hand, a strong ferromagnetic coupling 
between localized moments and itinerant electrons may 
arise in the transition metal systems due to the strong 
Hund’s coupling. In the strong coupling limit, the spins 
of itinerant electrons are fully polarized along the local¬ 
ized moments. Such a situation is well described by the 
double-exchange (DE) model. In this model, the kinetic 
motion of electrons induces an effective ferromagnetic 
(FM) interaction between localized momentsi^i^ This is 
called the DE interaction, which stabilizes a metallic 
FM state at low temperature. In these oxides, however, 
there generally exists antiferromagnetic (AFM) super¬ 
exchange (SE) interaction between the localized moments 
as well, and the competition of two interactions may give 
rise to nontrivial phenomena. One such example is found 
in perovskite manganese oxides. They are renown for the 
colossal magnetoresistance^^— in which inhomogeneity in 
the competing region between FM and AFM phases have 
been studied in relation to the magnetoresistancejiSrJi 

While the simple picture based on the above argument 
appears to work well in the manganese oxides, recent 
numerical studies of DE models on geometrically frus¬ 
trated latticesi^A^ have discovered emergence of inter¬ 



FIG. 1. (Color online). Schematic picture of a pyrochlore 
lattice. Each figure indicates the two-spin exchange interac¬ 
tions, Ji, and the four-spin interactions, Ki {I is an integer). 
For the latter, four Ising spins at the two colored sites (fc, 1) 
and two black sites {i,j) interact with each other. The arrows 
on the sites in the right panel indicate anisotropy axes of the 
Ising moments, rii. See the text for details. 


mediate phases in the competing region. In the case of 
checkerboard and triangular lattices, instabilities toward 
noncoplanar spin orderings were observed in the Monte 
Carlo (MC) simulation for the models with localized clas¬ 
sical Heisenberg momentsi ^^d^ Meanwhile, a thermally- 
induced intermediate phase with spontaneously-broken 
spatial inversion symmetry was found in the model with 
Ising moments on a pyrochlore latticeJ^ These results 
imply that, in the frustrated systems, the subdominant 
interactions beyond the simple DE mechanism poten¬ 
tially give rise to the nontrivial phases in the phase 
competing region. Indeed, in the previous study by the 
authors effective further-neighbor interactions derived 
from a strong coupling theory successfully predicted the 
presence of the intermediate phase. However, the method 
of the strong coupling expansion was not described in de¬ 
tail, and its accuracy has not been examined systemati¬ 
cally. 

In this paper, we present the framework of the strong 
coupling theory for the effective spin interactions between 
the localized moments in DE models. We illustrate the 
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technique for a spin-ice type DE model on a pyrochlore 
lattice (Fig. [T|). We show that the second-order expan¬ 
sion gives rise to various four-spin interactions in addition 
to the two-spin interactions between second- and third- 
neighbor sites. The results for the effective spin model 
are compared with numerical results for the DE model, to 
test the accuracy of this theory. We find that the expan¬ 
sion up to second order correctly reproduces the trend of 
magnetic phases while changing the AFM SE interaction 
and the itinerant electron density. 

The organization of this paper is as follows. In 
Sec. mi we introduce the model we consider, and elab¬ 
orate the details of strong coupling expansion used here. 
In Sec. uni we investigate the accuracy of this method 
by comparing the results with those of numerical diago- 
nalization and MC simulation. In the last, an extension 
of our theory to Heisenberg and XV localized moments 
is briefly discussed in Sec. IlYl Section El is devoted to 
discussions and summary of this paper. 


II. MODEL AND METHOD 

In this section, we introduce the model and method 
we used. In Sec. Ill A[ we introduce the DE model we 
consider in this paper. We explain the strong coupling 
expansion method in Sec. Ill Bl 


A. Model 


The DE model we consider in this paper consists of 
itinerant electrons and classical localized moments that 
are strongly coupled to each other. The Hamiltonian is 
given in the general form^ 

H = — -t- - JijSi ■ Sj, (1) 


where Ci (cj) is the annihilation (creation) operator of an 
itinerant electron and St is the localized moment at ith 
site. The electrons are described by spinless fermions, as 
their spins are perfectly polarized parallel to the local¬ 
ized spins. The first sum is the kinetic term of itinerant 
electrons. The transfer integral depends on the rel¬ 
ative position of ith and jth sites, and also depends on 
the localized spins at the two sites 


, _ .(bare) 

''d “ ‘'ij 


0i 0j 

COS — COS — 

2 2 

+ sin — sm exp[i((pj 



( 2 ) 


ith and jth sites. In this paper, we particularly consider 
the case of Ising local moments, i.e.. Si = with 
being a unit vector parallel to the anisotropy axis at ith 
site. 


B. Strong coupling expansion 


To investigate the magnetic properties of the model in 
Eq. ([T]), understanding the nature of effective spin inter¬ 
actions mediated by the coupled fermions is of crucial 
importance. To evaluate the effective spin interactions, 
we start by approximating the transfer integral in Eq. (ED 
by its amplitude. 
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where Oij is the angle between Si and Sji^ The DE 
model with this approximated form of the transfer inte¬ 
gral, tij, was used to study the magnetic properties of 
manganese oxides in several previous worksi^^— Using 
this approximation, and considering that the localized 
moments are of Ising type, we can rewrite the transfer 
integral as 

~ ( 4 ) 


where Si = Si ■ rii = ±1 is the projected spin parameter 
along rii, and 
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are real coefficients. Here, Ofj is the relative angle be¬ 
tween rii and rij. 

We consider the hopping term with the coefficient 
as the unperturbed Hamiltonian, Hq, and perform the 
expansion of Matsubara Green’s function with respect to 
the remaining term with tk. Hi. The Dyson equation is 
given by 


Gij{iuj) = gi,j{iu!) 


k,l 


3i,fe(iw) 


tlASi 


Gij{iuj). 


( 7 ) 


Here, Gij{iu!) is Matsubara Green’s function and gij(iui) 
is bare Green’s function of the unperturbed Hamiltonian. 
The term in the square bracket in Eq. ED is the scattering 
by Hi. The internal energy of the system is given by 


Here, 6i and (j)i are the polar and azimuthal angles of Si, 
respectively; is the transfer integral between ith 

and jth sites in the absence of the coupling to localized 
spins. The second sum in Eq. CD is the AFM SE interac¬ 
tion term in which is the exchange coupling between 


E— -l- ^'^^JijSi ■ Sj. (8) 

i,i i,i 

By replacing (c-Cj) by and expanding 

Green’s functions using the Dyson equation, one obtains 
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the effective spin model: the energy in Eq. ([8]) gives the 
effective Hamiltonian for the Ising spins Si. In this paper, 
we consider the expansion up to which leads 

to effective four-spin interactions in addition to two-spin 
ones. 

Regarding to the accuracy of this method, we note that 
in Eq. (0, \tS\ becomes small when is close to tt/2. 
Hence, it is expected that the perturbation is expected to 
be accurate when 0?^ ~ 7r/2, namely, the local anisotropy 
axes are perpendicular to each other. On the other hand, 
the approximation becomes less accurate as we approach 
the collinear case, = t},. 

In the following sections, we test this method for a DE 
model on a pyrochlore lattice with only nearest-neighbor 
(NN) transfer integrals and the localized Ising moments 
having spin-ice type anisotropy] ^^’^^ In this model, the 
anisotropy axes of two NN spins have the relative angle 
of ~ 109°, which is close to 7r/2. 


III. RESULTS 


In this section, as the benchmark of the method in the 
previous section, we study the effective spin interactions 
in a spin-ice type DE model on a pyrochlore lattice. In 
Sec. IIII A[ we present the effective spin model obtained 
from the strong coupling expansion. The accuracy of this 
method is investigated by comparing the ground state 
energy (Sec. IIII HI) and magnetic phase diagram for n 
and J fSec. IIII Cl) . The later is obtained by a variational 
method. The relevance of the phase diagram is further 
investigated in Sec. IIII Eil bv a MC method. 


A. Effective spin interactions 


To test the accuracy of the strong coupling theory, 
we consider an Ising spin DE model on a pyrochlore 
lattice with hoppings and SE interactions for NN sites 
only. The anisotropy axes for Ising moments are de¬ 
fined along local [111] axes: ria = (l/-\/3)(l, —1, — 1), 
rih = (l/y3)(-l, 1,-1), ric = (1/V5)(-1,-1,1), and 
«d = (1/V3)(1,1,1) (see Fig. d]). For simplicity, we re¬ 
place the NN hopping integral by iij as in Eq. O; the 
starting Hamiltonian is given by 

H = - + H.C.) + J ^ 5. • S„ (9) 

(hi) (iJ) 


where the sum (i,j) is taken over the NN sites on the 
pyrochlore lattice. Hereafter, we take the bare hopping 
for NN sites, t = 1, as the energy unit. 

Applying the strong coupling expansion in Sec. Ill Bl to 
this model, we construct an effective Ising model on the 
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FIG. 2. (Color online). Effective spin interactions mediated 
by itinerant electrons: (a) two-spin and (b) four-spin interac¬ 
tions as functions of the electron density n. The definition of 
interactions is given in Fig. [T] 


pyrochlore lattice, whose Hamiltonian is given by 

hj {i,j),{k,l), 

( 10 ) 

up to a constant given by the contribution from Hq. Here, 
the first term represents effective two-spin interactions 
for nearest-, second- and third-neighbor sites, while the 
second term describes the four-spin interactions. The 
last term is the AFM SE interaction already present in 
Eq. (|ni), where the coefficient —1/3 comes from the pro¬ 
jection of the Ising spins to the anisotropy axes at each 
site. In the model in Eq. ®, the leading order in ex¬ 
pansion gives NN two-spin interaction while the second 
order gives four spin interactions between spins on the 
edge of two bonds at arbitrary distance. The second and 
third neighbor two-spin interactions also arise from the 
second order expansion, when the two bonds share a site, 
e.g., j = k. In general, the two-spin interactions between 
sites with / Manhattan distance arise from the Ith. order 
in expansion, i.e., they decay exponentially with distance 
in contrast to the RKKY interaction. 

Figure H] shows the effective spin interactions that arise 
from the expansion for the model in Eq. m as func- 
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FIG. 3. (Color online), n dependence of the ground state 
energy for different spin configurations at J = 0. The lines are 
the results obtained by numerical diagonalization of the DE 
model in Eq. and the symbols are those for the effective 
spin model in Eq. uni obtained from the strong coupling 
expansion. See the text for details. 


FIG. 4. (Color online). J dependence of the ground state en¬ 
ergy for different spin configurations calculated for the effec¬ 
tive spin model in Eq. GQl) obtained from the strong coupling 
theory [(a) and (b)] and the DE model in Eq. ([51 by numer¬ 
ical diagonalization [(c) and (d)]. The energies are measured 
from those for the 2120 state. The electron density were set 
at n = 0.3 (0.4) for (a) and (c) [(b) and (d)[. 


tions of the itinerant electron density n = '^i{c\ci)/N 
{N is the number of sites). We here show the results 
for 1/4 < n < 1/2, as the higher order terms appear 
to be relevant for lower n (see below) and the flat band 
may lead to some complexity for higher n. The defini¬ 
tion of each interaction is shown in Fig. [T] As shown in 
Fig. [21 the most dominant interaction mediated by itin¬ 
erant fermions is the FM NN interaction Ji, consistent 
with what is expected in DE models. This dominant in¬ 
teraction gives rise to instability toward a FM state when 
the SE interaction J between localized moments is suffi¬ 
ciently weak. 

In addition, we also have other two-spin and four-spin 
interactions which are an order of magnitude smaller 
than Ji. These interactions potentially become impor¬ 
tant in the phase competing region between DE-driven 
EM and SE driven AFM phases, i.e., when the SE in¬ 
teraction cancels the NN FM interaction. Some of these 
subdominant interactions are also plotted in Fig. [2] In 
this density range, most of the four-spin interactions de¬ 
cay rapidly with distance; Kq is the interaction of four 
spins on the same tetrahedron, while (i ^ 0 ) are inter¬ 
actions between spins with further distance. As a conse¬ 
quence, the dominant interactions are J 2 , J 3 , and iFo for 
electron densities n > 1/4. The result implies that, for 
n > 1/4, only considering a limited number of interac¬ 
tions are sufficient in reproducing the qualitative nature 
of the model. 


ders with all spins on a tetrahedron pointing inward or 
outward (all-in/all-out; AIAO) and two spins inward and 
two spins outward (two-in two-out; 2120 ) as an example 
of the magnetic order with AFM and FM NN correlation, 
respectively. In addition to these phases, two magnetic 
orders that were found in the related Kondo lattice mod¬ 
els, 32-sublattice (32-sub) order— and spin cluster (SC) 
orderare also considered. 

Figure El shows the result of the ground state energy 
calculated for different magnetic orders in the effective 
spin model in Eq. (fTUl) by taking into account Ji, J 2 , 
J 3 , and Kq. Eor n > 1/4, the results of strong coupling 
expansion are in accordance with those of numerical di¬ 
agonalization for the model in Eq. (|9]). As the NN EM 
interaction Ji gives the largest contribution to the ground 
state energy, this result shows that the estimate of Ji is 
successful in the strong coupling expansion. 

We also calculated the ground state energy for n < 
1/4. In this region, however, the results show strong 
deviation from the result of numerical diagonalization 
(not shown here). This is presumably due to the pres¬ 
ence of longer-range interactions than second- and third- 
neighbors which are more dominant than the case for 
n > 1/4. 


C. Variational phase diagram 


B. Numerical diagonalization 

To examine the accuracy of above theory, we start by 
evaluating the ground state energy. The results are com¬ 
pared to those of numerical diagonalization for several 
different spin configurations. We here consider q = 0 or- 


Next, to evaluate the effect from further-neighbor two- 
spin interactions and four-spin interactions in Eq. m, 
we study the ground state phase diagram by a variational 
calculation with changing the AEM SE interaction J. As 
J tends to cancel the effective NN FM interaction from 
the DE mechanism, this essentially corresponds to study¬ 
ing the effect of subdominant interactions that arise from 
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FIG. 5. (Color online). Variational phase diagrams for (a) the effective spin model obtained by the strong coupling expansion 
[Eq. (llUll l. (b) the DE model in Eq. (O, and (c) the DE model in Eq. ([T]). The phase diagram in (a) is calculated by only 
considering Ji, J2, J3, and Kq. 


the second order expansion in the strong coupling theory. 
For this purpose, we here perform a variational calcula¬ 
tion comparing the ground state energy for different spin 
configurations: 2120, 32-sub, SC, and AIAO orders^i 

Figure m shows the results of ground state energy for 
different magnetic orders measured from that for the 
2120 state. Figure SKa) shows the result for the effective 
spin model in Eq. (HOD obtained by the strong coupling 
theory at n = 0.3. We find all four phases with increas¬ 
ing J; the transition takes place from the 2120 to 32-sub 
state at J ^ 0.109, to SC state at J ~ 0.121, and to 
AIAO state at J ~ 0.157. A similar trend is also found 
in the result calculated by numerical diagonalization of 
the model in Eq. (0) [Fig- Sic)]. Note that the electronic 
phase separation is not considered here; see Fig. [S] 

In the results for n = 0.4, the strong coupling theory 
predicts four magnetic phases in the ground state, as 
shown in Fig. IHb). On the other hand, as shown in 
Fig. lUJd), there are only three phases for the model in 
Eq. (0): 2120, 32-sub, and AIAO phases. 

Performing the above variational calculation while 
varying n and J, we map out the ground state phase dia¬ 
gram. Figure0shows the variational phase diagrams cal¬ 
culated for the effective spin model and the DE models. 
Figure IHKa) shows the result for the effective spin model 
constructed from the strong coupling theory [Eq. (fTOl) ]. 
For all 1/4 < n < 1/2, we found both SC and 32-sub 
phases between the 2120 and AIAO phases. 

On the other hand, the DE model shows a slightly 
different phase diagram. Figure |5][b) shows the phase 
diagram for the model in Eq. (0). In Fig. [51[b), we found 
only the 32-sub state in the intermediate range of J for 
0.321 < n < 0.466, while only the SC state appears for 
n < 0.288 and n > 0.471; the two phases appear in the 
narrow range of 0.288 ^ n < 0.321 successively while 
increasing J. Between these phases as well as to the 
2120 and AIAO phases, the system exhibits an electronic 
phase separation, which is due to the discontinuity in n 
associated with the first-order magnetic phase transition. 
The result indicates that, with varying n, there are two 
regions in the phase diagram where the phase competing 


region is dominated by either 32-sub or SC order. This 
result is in contrast to the strong coupling theory, where 
the two phases appear for all n. Nevertheless, the strong 
coupling theory predicts the trend of instability toward 
the correct intermediate ground states found in the DE 
model. 

In the last, we discuss the variational phase diagram for 
the DE model in Eq. o, without approximating by tij 
in Eq. O- The result is shown in Fig.[Sl[c). The SC phase 
appears only for n < 0.25, while the 32-sub phase appears 
for 0.323 < n < 0.418; no intermediate phase is found for 
0.25 ^ n < 0.323 and n > 0.418. Hence, despite the 
absence of quantum phase in the hopping term, the two 
DE models show qualitatively similar ground state phase 
diagrams. This is a crucial observation for justifying the 
approximation we used in this paper, as we ignored the 
effect of the quantum phase in the strong coupling theory. 

D. Monte Carlo simulation 

For further comparison, we study the models in Eq. (ED 
using the unbiased MC method, which allows calculat¬ 
ing thermodynamic quantities of the model in Eq. ED 
without any approximations)^ This method and its vari- 
ants^Si^l have recently been used to explore unconven¬ 
tional phases in the phase competing region of the DE 
models on frustrated lattices.— “— 

The calculations were done with the system size N = 
4 X Ns with Ns = 4^ under the periodic boundary con¬ 
ditions. Thermal averages of physical quantities were 
calculated for typically 3600 MC steps after 600 steps for 
thermalization. Some of the low-temperature data were 
calculated for longer MC steps up to 10400 steps. We di¬ 
vided the MC measurements into five bins and estimated 
the statistical errors by the standard deviations among 
the bins. 

MC results for both DE and effective spin models at 
n = 0.25 have already been published by the authors^ 
which are in accordance with the variational phase dia¬ 
gram in Fig. EKc). In the DE model, only the SC state 
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FIG. 6 . (Color online). Results of Monte Carlo simulation for 
the DE model in Eq. © at n ~ 0.38: (a) p 22 , (b) p 40 , and 
(c) S{q)/N. See the text for details. 


was found at the lowest temperature while increasing J, 
as the ground state in intermediate region. In addition, 
a fluctuating SC state with spatial inversion symmetry 
breaking was obtained in the intermediate temperature 
range above the critical temperature of the SC order. 
However, we did not hnd the 32-sub order in accordance 
with the variational phase diagram. 

To check the existence of the 32-sub order in the region 
n > 1/4, we here performed the MC simulation for the 
DE model in Eq. (©). Figure [ 6 ] shows the result of MC 
simulation for iV = 4 x 4^ sites at n = 0.38. Figures (HKa) 
and llKb) show ratio of tetrahedra with two-in two-out 
{ 922 ) and all-in/all-out (p 4 o) configurations. When J = 
0.06, the results show enhancement of 922 and decrease 
of P 40 with decreasing temperature. This is a sign of 
the EM correlation for the NN bonds i^^i^^ However, we 
could not reach the transition temperature for the long- 
range order in our simulation due to the freezing at lower 
temperatures. In contrast, for J = 0.22, p 4 o increases 
with decreasing temperature, approaching p 4 o —?► 1. In 
addition, the structure factor 

S'a(q) = ;^ XI 

i^jGcx 

for q = 0 shows increase from zero as shown in Fig. Elc), 


indicating the AIAO order in the ground state. 

On the other hand, the result for J = 0.14 shows 
922 —f 0.375 and /O 40 —>■ 0.125 with decreasing temper¬ 
ature [Figs. EKa) andEKb)]. This is a sign of the phase 
transition to 32-sub order, where 6/16 (2/16) tetrahe¬ 
dra are in 2120 (AIAO) spin configuration!^ This tran¬ 
sition is also confirmed by the spin structure factor for 
q = ( 7 r, 7 r, 7 r) plotted in Fig. EKc). 

From these results, we conclude that, a sequence of 
states, 2120, 32-sub, and AIAO, appears for n = 0.38 
with increasing J. On the other hand, at n = 0.25, the 
presence of SC order in the intermediate J region was 
previously reported)^ consistently with the variational 
phase diagram in Fig.EKc). The results support that the 
effective spin model by the strong coupling theory can 
predict the trend of intermediate phases in the competing 
region in the DE model. 


IV. EXTENSION TO HEISENBERG AND XY 
MOMENTS 


In the last, we briefly discuss potential extension of 
this method to DE models with Heisenberg or XY type 
localized moments. In these cases, we cannot use the 
expansion in Eq. (©). One alternative approach is to 
replace cos 0ij in Eq. (©) by • Sj, and consider — 
^^(bare) perturbation, where 0 < a < 1. In the 

leading order, this gives EM NN interaction of the form 
\Jij\ oc tij] the classical spin model with Jij as the NN 
interaction was studied motivated by the DE modelsi^ 

Another route to evaluate the effective interactions is 
to use an asymptotic expansion, for instance, Taylor se¬ 
ries 




S^ ■ Sj, 




S, ■ S, {S, ■ 


( 12 ) 

(13) 


This expansion naturally predicts the presence of posi¬ 
tive biquadratic interaction as the subleading interaction. 
Recently, the positive biquadratic interaction has been 
proposed^^ as the possible origin of non-coplanar phases 
found in the weak coupling region of a triangular Hondo 
lattice modeli^^— The same non-coplanar state has also 
been found in the DE limit, but its mechanism have not 
been explained so farJ^ A simple argument based on our 
theory suggests that the non-coplanar phase in the DE 
limit may also be stabilized by the biquadratic interac¬ 
tion that arises from the strong coupling expansion. This 
also implies that, in the DE models with Heisenberg mo¬ 
ments, even 0 (tk) in the expansion may give rise to un¬ 
conventional magnetism. 
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V. DISCUSSIONS AND SUMMARY 

To summarize, in this paper, we studied the fermion- 
mediated effective spin interaction in an Ising spin 
double-exchange model on a pyrochlore lattice. To eval¬ 
uate the effective interactions from microscopic theory, 
we used a strong coupling expansion and calculated the 
effective interactions up to second order in terms of the 
spin-dependent electron hopping. We showed that, ef¬ 
fective four-spin interactions appear in the second order 
expansion. We also found that the effective interactions 
are limited to short range for this model, at least, for 
electron density 1/4 < n < 1/2. 

Focusing on this region, we studied the accuracy of 
strong coupling theory. From comparison to the nu¬ 
merical results on the double-exchange model, we found 
the strong coupling theory gives a good estimate of the 
ground state energy. In addition, we studied the ground 
state phase diagram in the presence of antiferromag¬ 
netic super-exchange interactions between the localized 
moments. The calculations were done by a variational 
method, and some of the results were also confirmed by 
a Monte Carlo simulation. We found that the strong cou¬ 
pling method correctly captures the trend of intermediate 
phases found in the double-exchange models. 

It is interesting that the expansion up to second order 
appears to capture the correct trend of the DE model. 


although the energy scale we are discussing is very small, 
E/t ^ 0.01. This may be related to the fact that the 
ground state of a spin model is often sensitive only to the 
sign and not to the magnitude of interactions. For the ef¬ 
fective spin model we studied, the 32-sublattice and spin- 
cluster states appear in the region of both second- and 
third-neighbor exchange interactions being antiferromag- 
neticJ^ Hence, although the second order in expansion 
is not sufficient in correctly predicting which phase wins, 
it is still successful in guessing the correct candidates for 
the intermediate phase. This arguments cast a possi¬ 
ble restriction on the application of the strong coupling 
theory we proposed; the strong coupling theory predicts 
correct trend only when the double-exchange model has 
a few subdominant interactions that are relevant, in ad¬ 
dition to the ferromagnetic nearest-neighbor interaction. 
However, if this condition is satisfied, we can expect that 
the theory gives the correct trend. 
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